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Effects of the density dependence of the nuclear symmetry energy on ground-state properties of 
superheavy nuclei are studied in the relativistic mean-field theory. It is found that the softening 
of the symmetry energy plays an important role in the empirical shift [Phys. Rev. C 67, 024309 
(2003)] of spherical orbitals in superheavy nuclei. The calculation based on the relativistic mean- 
field models NL3 and FSUGold supports the double shell closure in ^^^120 with the softening of 
the symmetry energy. In addition, the significant effect of the density dependence of the symmetry 
energy on the neutron skin thickness in superheavy nuclei are investigated. 
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I. INTRODUCTION 

The persistent interest in the synthesis of superheavy nuclei (SHN) acquires refreshment due to 
recent progresses [1-14]. This is a hot field where people often expect the next new superheavy 
element (SHE) that can be synthesized in the laboratory. Indeed, since the cross section of the SHE 
synthesis is very small, it is much more difficult to synthesize the heavier and heavier SHE [2] . For 
instance, the cross section of the cold fusion reduces almost exponentially with the increase of the 
nuclear charge in the superheavy region. One of important factors that affects the synthesis is the 
shell closure in superheavy nuclei. However, predictions turned out to be quite divisive for various 
theoretical approaches. For instance, the microscopic-macroscopic model predicts the the double 
shell closure at (Z=114,N=184) [15]; various nonrelativistic models with Skyrme forces can predict 
different double shell closures at (Z=114,N=184) [16], (120,172) [16-18] or (126,184) [16, 17, 19], 
while most relativistic mean-field (RMF) models incline the double shell closure at (120,172) [16, 
17, 20]. In general, the diversity of predictions on the shell closure in the superheavy region is 
associated with various single-particle properties near the Fermi surface. 

In the recent decade, the extraction of the constraint on the density dependence of the sym- 
metry energy has been another hot spot in nuclear physics due to the availability of high-quality 
radioactive beams. The density dependence of the symmetry energy plays important roles in un- 
derstanding many important issues in astrophysics, see, e.g., Refs. [21-23], properties of proton- 
or neutron-rich nuclei and the reaction dynamics of heavy-ion collisions, see, e.g., Refs. [24-27]. 
However, the density dependence of the symmetry energy is still poorly known especially at high 
densities [27]. Recent extraction of the neutron skin thickness of ^''^Pb from collective flow data 
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of heavy ion collisions [28-30] exhibited the softening tendency of the symmetry energy. Since the 
density dependence of the symmetry energy can reflect the surface property of the isovector po- 
tential, the effect on the single-particle property and characteristic of the shell closure in SHN may 
be induced by the softening of the symmetry energy. Moreover, it was found that the existence of 
the central depression is important for the double shell closure in ^^^120 [16, 31, 32]. In presence 
of the central depression, the sensitivity of the properties of SHN to various density dependences 
of the symmetry energy can be afi'ected. Though the properties of SHN have been explored in a 
great number of works [15-20, 31-45], the investigation on the symmetry energy dependent effect 
is scarce. Thus, it is the aim of this work to investigate the effect of the softening of the symmetry 
energy on ground-state properties of SHN, especially the shell closure. 

In the past, the isoscalar-isovcctor coupling was first introdiiccd in RMF models to mimic various 
density dependences of the symmetry energy in Ref. [22], and its effects on the properties of finite 
nuclei, nuclear matter and neutron stars have extensively been investigated in the literature [22, 
24, 25, 46-50]. It is an economic way to simulate various density dependences of the symmetry 
energy with the inclusion of the isoscalar-isovector coupling in RMF models. In addition, the 
RMF theory is successful in describing the properties of nuclei from proton drip line to neutron 
drip line, since it can provide a dynamic description for the spin-orbit interaction, e.g., see reviews 
in Rcfs. [51 53]. In this work, we thus perform the investigation with RMF models. The paper is 
arranged in the following. In section II. the brief formulas are given for RMF models. The results 
and discussions are presented in section III. At last, a summary is given in section IV. 



where ^/i, a, w, and bo are the fields of the nucleon, scalar, vector, and neutral isovector-vector, with 
their masses Mj^,mcr,mu, and rup, respectively. A^j, is the photon field. gi{i = a,LO,p) are the 
corresponding meson- nucleon couplings. Fp^i,-, and A^^, are the strength tensors of w and p 
mesons, and photon, respectively 



II. A BRIEF FORMALISM 



The relativistic lagrangian can be written as: 



£ = ^[n^a'* -Mn + g^a - g^-y^uj" - gp^^nbi^ - e^(l + r3)7^^'']V 



(1) 



(2) 



The self-interacting terms of tr, w mesons and the isoscalar-isovector couphng are given generally 
as 

+AglglKw^u^'boX- (3) 

Here, the isoscalar-isovector coupling term is introduced to modify the density dependence of the 
symmetry energy. 

Using the Euler-Lagrangian equation, the equations of motion for nucleons and mesons can be 
obtained. In the RMF approximation, the mesons are approximated by their classic fields with 
quantum motion neglected. The Dirac equation in RMF is written as 

[-ia ■ V + /3Mj^ + g^uJo{r) + QpT^boir) +e^{l + r3)Ao(r)]Va(r) = E^Va, (4) 

with Mjljr = Mjv — 9(T(^{r) and the single-particle energy. For simplicity, the isospin subscript 
for the p-meson field is omitted hereafter. For the mesons and photon, the equations of motion are 
given as 

(A - m|)</)(r) = -s^r) (5) 

where for the photon, = 0, and 

9aPs{r) - g2C7^{r) - gacr^ir), a , 



(r) = 



(6) 



9wpB{r) - caw^ - d>glglKuJo{r)hl{r), w , 
ffpP3(r) - ^glglKho{r)ujl{r), rho , 

epc{r), photon. 
Here ps, Pb, Pa and pc are the scalar, vector, isovector and charge densities, respectively. We see 
that the fields bo and uq can be modified by the isoscalar-isovector coupling. This modification 
can also affect the spin-orbit potential which is written as 

where 

= M^f + E^- V{r) « 2Miv - (g^r) + V{r)), 
V{r) = gu;Mr)+9pbo{r)t + e{^-^)Ao, (8) 

with t = ±1 for the proton and neutron, respectively. 



The total binding energy is given as 



B = En + + -Bt^o + Eba + Ec + Ecm 

a 



+ 2, 

+ \9P j d^rho{r)[p:i{r)+SgpglKi^l{r)bo{r)] 

+ \e j dVAo(r)pe(r) - ^41^1/3. (9) 

In practical calculations, the BCS pairing interaction is also included using the constant pairing 
gaps which arc obtained from the prescription of Moller and Nix [54]: A„ = 4.8/-/V-'^/^, Ap = 
4.8/,Z^/^ with N and Z the neutron and proton numbers, respectively. This prescription was also 
used for the SHN in Ref. [43]. The cut-off 82^4"^/^ MeV above the nucleon chemical potentials is 
used to normalize the pairing energy [52]. The coupled Dirac and meson equations are solved for 
spherical nuclei with an iterative procedure. Details in solving the equations can be found in the 
literature [51-53], and are not reiterated here. 

III. RESULTS AND DISCUSSIONS 

We first study the properties of the SHN with the RMF parameter set NL3 [55] where the 
isoscalar-isovector coupling is taken into account to mimic various density dependences of the 
symmetry energy. For comparisons, calculations are also performed with the RMF parameter set 
FSUGold [48] that features the isoscalar-isovector coupling. In RMF models, the symmetry energy 
can be written as: 

I ( gp\^ kl 1 , fc| 



where m* is the p-meson effective mass with m* = -|- d>K^{gugpUJoY, 5 is the isospin asym- 
metry with S = Ps/pB, and Ep is the Fermi energy. The first term is the potential part of the 
symmetry energy, and the second term is the kinetic part. The modification to the symmetry 
energy is dictated by the potential part through the isoscalar-isovector coupling. For a given Av, 
we follow Ref. [22] to readjust the pNN coupling constant gp so as to keep the symmetry energy 
unchanged at kp = 1.15 fm~^ (p = 0.7po)- As shown in Fig. 1, the symmetry energy is softened 
by the isoscalar-isovector coupling. With this softening of the symmetry energy, the appreciable 
reduction of the neutron skin thickness in heavy nuclei can be obtained without compromising 
the success in reproducing a variety of ground-state properties [22] . Due to the inclusion of the 
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isoscalar-isovector coupling, the total binding energy of heavy nuclei changes by a few MeV, and in 
SHN this change can rise moderately. To reduce the variation of the binding energy in SHN, one 
may readjust slightly the parameters such as the meson-nucleon coupling constants and mesons. 
Without priority, here we readjust slightly the a meson mass m.rr. For simplicity, we do not perform 
the best-fit procedure, and the value of is just refitted to the binding energy of ^°^Pb. The 
readjusted parameters with various Ay and properties of ^°^Pb are listed in table I. Except for 
the original parameter sets NL3 and FSUGold, other parameter sets listed in Table I are named 
according to the value of Ay. Next, we perform calculations for SHN with these parameter sets and 
examine the sensitivity of ground-state properties of SHN to differences in the symmetry energy. 
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FIG. 1: Density dependence of the symmetry energy with various isoscalar-isovector couplings in NL3 and 
FSUGold. 

In Fig. 2, wc plot the single-particle energies for ^^^120 in the NL3 calculations. Results with 
various isoscalar-isovector couplings are displayed in columns. The large gaps for N=172 and 
Z=120, as shown in Fig. 2 indicate that the nucleus ^^^120 is doubly magic. It is seen that the shell 
closure at N=172 and Z=120 undergoes a small but favorable enhancement due to the inclusion of 
the isoscalar-isovector coupling. As shown in the left panel, the position of 7rl/i9/2 relative to that 
of 7r3si/2 shifts appreciably with the inclusion of the isoscalar-isovector coupling. For the large 
Ay, even the level inversion takes place. This shift is in favorably agreement with the one called 
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TABLE I: Readjusted parameters in NL3 and FSUGold with ground-state properties of Pb. The binding 
energy per nucleon (B/A), proton radius (vp) and neutron skin thickness {rp — r„) are listed. The shghtly 
modified incompressibihty is listed in the last column. 



Model 


Av 


9p 


m„ (MeV) B/A (MeV) 


Vp (fm) 


r-n - rp (fm) 


K (MeV) 


NL3 


0.000 


4.4740 


508.194 


7.889 


5.459 


0.281 


271.78 


NL3wl5 


0.015 


4.9652 


508.240 


7.890 


5.465 


0.238 


272.25 


NLSwSO 


0.030 


5.6642 


508.270 


7.890 


5.475 


0.195 


272.56 


NL3w50 


0.050 7.3236 


508.270 


7.890 


5.496 


0.132 


272.56 


FSUGwl5 


0.015 


5.0403 


491.490 


7.883 


5.463 


0.248 


229.96 


FSUGold 


0.030 


5.8837 


491.500 


7.883 


5.473 


0.207 


230.00 


FSUGw45 


0.045 


7.3695 


491.480 


7.883 


5.488 


0.158 


229.92 
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FIG. 2: Single-particle energies in 120 with various parameter sets in the NL3 calculations. 

as the empirical shift in Refs. [33, 34]. Up to now, there is no direct data of the single-particle 
energies of SHN, while the so-called empirical shift is obtained by extracting available single-particle 
energies in deformed nuclei of the A ^ 250 region (e.g., ^^^Bk) [56-59]. Since several deformed 
single-particle levels observed in the A ~ 250 nuclei emerge from spherical subshells in SHN, the 
appropriate description of empirical shifts can provide a favorable support for the predictions on 
properties of SHN, especially the nuclear magicity. With the inclusion of the isoscalar-isovector 
coupling, the empirical shift between the 'jTlhg/2 and 7r3si/2 can be well reproduced. It is also 
interesting to see that the low-j levels 7r3p3/2 and 7r3pi/2 can be significantly modified by the 
isoscalar-isovector coupling. However, its influence on the shell closure at Z=120 remains small. 
As a result, the shell gap for Z=120 is just weakly affected by the empirical shift. This means 
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that these parameter sets including the original NL3 can provide a reliable prediction on the shell 
closure at Z=120. On the other hand, the empirical shift for [33] in the single-neutron 

spectrum is not reproduced with the inclusion of the isoscalar-isovector coupling. In Ref. [33] , it 
can be scon that the shell closure at N=172 is just moderately affected by the empirical shift. The 
present prediction on the large gap for N=172 does not contradict with the early analysis with 
the empirical shift. Indeed, the relative energies between v2ggi2, v^g-jji and 1/3^5/2, one of which 
determines the N=172 gap, are almost independent of RMF parametrizations, see Ref. [33] and 
references therein. The situation of the neutron shell closure at N=184 is less known, since there 
is no empirical constraint on the i^4si/2- However, the N=184 gap is almost independent of the 
shifts of the interior levels caused by the inclusion of the isoscalar-isovector coupling. Similarly, 
the N=184 gap would not be much affected even if the empirical shift for is accurately 

reproduced. In this sense, the occurrence of the shell closure at N=184 seems unlikely. 



proton neutron 




Radius (fm) 

FIG. 3: Nucleoli potentials (upper panels) and nucleon density distributions (lower panels) in '^^^120 with 
various parameter sets in the NL3 calculations. The nucleon potential is defined as ?7 = V{r) — ga<j{r), 
also see Eq.(8). 

Now, let's understand underlying factors that intrigue the significant shift in single-particle spec- 
tra with the inclusion of the isoscalar-isovector coupling. As shown in Fig. 2, the spin-orbit splitting 
can be modified by the isoscalar-isovector coupling, and the modification increases moderately with 



the angular momentum. However, the modification to the spin-orbit coupling is just moderate and 
is not sufficient to cause the empirical shift. Another factor that affects the modification to the 
single-particle spectrum is the orbit-orbit interaction. The orbit-orbit interaction can generally 
be given in a form of the centrifugal force, and it reflects the flatness of the nuclear potential. 
In Fig. 3, the nucleon potentials and density distributions in ^^^120 are plotted. As shown in 
the upper left panel of Fig. 3, the homogeneity of the proton potential in the central region can 
be modified by the isoscalar-isovector coupling. This modification can bring about the moderate 
change in the single-proton levels. Both changes in the orbit-orbit and spin-orbit interactions thus 
lead to significant empirical shifts of proton levels, observed in ^^^120. Similarly, the shifts in the 
single-neutron levels can be understood by the modification in the orbit-orbit and spin-orbit in- 
teractions caused by the isoscalar-isovector coupling. Though the isoscalar-isovector coupling can 
affect shifts in the single-particle levels for both protons and neutrons significantly,its modification 
to the proton or charge radius is much less than that to the neutron radius. This can be seen in 
Fig. 3 by comparing the modification to the neutron density distribution with that to the proton 
one. 

Furthermore, it is significant to establish a quantitative correlation between the density de- 
pendence of the symmetry energy and the relative shift of single-particle energies. It is known 
that the symmetry energy can be expanded in the vicinity of saturation density in the following 
form [60, 61]: 



J^symKPB) — J^symKPo) + 1 TT, 2 ' ' \^^) 

6 po ^o pQ 



where L and Ksym are the slope and curvature of the symmetry energy at saturation density, 
respectively, defined as 

dEs 



L = 3po ^ 



dpE 



n2 p 
2 ^sym 



9po « o 



(12) 



The slope of the symmetry energy defines the symmetry pressure through the relation Psym = 
poLjZ. In the history, the proper inclusion of the spin-orbit coupling played an important role in 
giving rise to the correct shape of nuclear potential and hence the ordering of the energy levels. 
Since the spin-orbit interaction is associated with the surface property subject to the subsaturation 
density region, it is useful to similarly define the slope and curvature of the symmetry energy at half 
saturation density, {Psym) and Kg^^. As seen in Fig. 2, the relative shifts of the orbitals change 
with respect to the isoscalar-isovector coupling. Since the correlation between these relative shifts 
and the density dependence of the symmetry energy are similar, we just plot as an example in Fig. 4 
the relative shift between 7rl/i9/2 and 7r3si/2 as a function of the symmetry pressure and curvature. 



9 



As shown in Fig. 4, the relative shift of these two levels is approximately linear in the symmetry 
pressure (at saturation density) and correlates quadratically with the symmetry pressure at half 
saturation density. A stronger correlation at half saturation density reflects the strong dependence 
of singlc-particlc energies on the surface property of finite nuclei. In the right panel, it is shown 
that the relative shift of the two levels is linear in the curvature at half saturation density, which 
is consistent with the quadric correlation as shown in the middle panel. It was found in the early 
time for the RMF theory that the proper density dependence of the potential is important for a 
correct spin-orbit potential and hence the ordering of orbitals, e.g. see [62, 63]. Similarly, it is 
here interesting to see that the density dependence of the symmetry energy (or of the isovector 
potential) affects moderately the single-particle energies. In particular, in the present work we 
can obtain the relative shift between 7rl/ig/2 and 7r3si/2 with the appropriate density dependence 
of the symmetry energy. The correlations shown in Fig. 4 also exist for other SHN, even if the 
central depression disappears. Here, the central depression plays a role in enhancing the correlation 
strength ( e.g. the slope in the linear correlation). 




2.5 5 1 1.5 -100 

PSYM (Me^/f"^"^) I^YM Ks^M(MeV) 

FIG. 4: Relative shift between TTlhg/2 and 7r3si/2 8s a function of the symmetry pressure and curvature. 
The line in the left and right panels are obtained with a linear fit. The pressure in the left panel is evaluated 
at saturation density, while the pressure and curvature in the middle and right panels, respectively, are 
calculated at half saturation density. 

Next, we turn to the discussion on the two-nucleon gaps. Besides the gap in the single-particle 
spectrum, a direct measure of the shell closure is the appearance of the peak in the two-nucleon 
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FIG. 5: Binding energies per nucleon of Z=120 isotopes (upper panel) and two-neutron gaps (lower panel) 
with various parameter sets in the NL3 calculations. 

gaps, which are defined as [17]: 

52n = S2n{N + 2,Z)-S2n{N,Z) = 2B{N,Z)-B{N-2,Z)-B{N + 2,Z), (13) 
S2p = S2piN,Z + 2)-S2p{N,Z) = 2B{N,Z)-B{N,Z-2)-B{N,Z + 2), (14) 

where S2n and S2p are the two-neutron and two-proton separation energies, respectively. The peak 
of the two-nucleon gap reflects the large change of the two-nucleon separation energy, signaling 
the shell closure. Moreover, the two-nucleon gap can reflect appropriately the gap size in the 
single-particle spectrum [16]. In Fig. 5, we plot the binding energy per nucleon (upper panel) and 
the two-neutron gap (lower panel) for the Z=120 isotopes in the NL3 calculations. As shown in 
the upper panel of Fig. 5, the difference between binding energies is small for different isoscalar- 
isovector couplings, and this is attributed to the refitting of rria- As shown in the lower panel of 
Fig. 5, the two-neutron gap at N=172 can earn a moderate rise with the inclusion of the isoscalar- 
isovector coupling, which is consistent with the observation of the single-neutron spectrum in Fig. 2. 
The similar increasing tendency also occurs for the N=198 shell gap. As shown in the lower panel 
of Fig. 5, the peak of the two-neutron gap occurs at N=172, 184 and 198 in the Z=120 isotopes. 
Comparing to the sharp peak at N=172, peaks at N=184 and 198 becomes much blunt. This 
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indicates that the shell gaps at N=184 and 198 are not well developed. Also, we examine the 62n 
for the Z=126 isotopes, and the case is similar. 

As shown in Fig. 5, the two-neutron gap can be affected by the isoscalar-isovector coupHng. In 
some isotopes such as ^^^120 and ^^*120, the two-neutron gap can gain a rise with the inclusion 
of the isoscalar-isovector coupling. However, the modification to the two-neutron gap caused by 
the isoscalar-isovector coupling is not governed by the isotopic effect in SHN. For instance, this is 
clear as comparing the modification to ^271 at N=172 with the one at N=184. In fact, the relatively 
pronounced modification is associated with the specific geometries such as the central depression or 
enhancement. It was pointed out in Ref. [32] that the central depression in ^^^120 is predominantly 
from the proton occupation of high-j orbitals, while the neutron central depression results from 
the strong coupling between protons and neutrons. In presence of the central depression, the 
change in the neutron potential and density distribution caused by the isoscalar-isovector coupling 
exhibits a radial inhomogeneity, as shown in right panels of Fig. 3. This inhomogeneity causes 
consistently modifications to the level shifts and the two-neutron gaps. With the increase of the 
neutron numbers, the neutron central depression tends to disappear, while the central enhancement 
appears with more neutrons. The inhomogeneity of the modifications produced in presence of the 
central enhancement explains the change in the two-neutron gap in more neutron-rich isotopes, as 
shown in the lower panel of Fig. 5. 
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FIG. 6: Two-proton gap 52p for N=172 isotones with various parameter sets in the NL3 calculations. 



To examine the proton shell closure, we plot in Fig. 6 the two-proton gap for the N=172 isotones. 
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As shown in Fig. 6, there is only one peak at Z=120. No peak is observed at Z=114 and 126. The 
proton shell closure at Z=126 was predicted by the Hartree-Fock approach with a variaty of Skyrme 
interactions [16, 17], while the RMF models disincline the appearance of this shell closure. As seen 
in Fig. 2, though the isoscalar-isovector coupling tends to shift the relative position between 7rliii/2 
and 7r3p orbitals, the formation of the Z=126 shell closure does not appear. In RMF models, the 
Z=114 shell closure was only predicted by parameter sets NL-SH and NL-RAl due to the relatively 
large spin-orbit splitting of 2/ orbitals [64, 65]. In NL3, the spin-orbit splitting of 2/ orbitals is not 
large, and though the modification of the isoscalar-isovector coupling to 52p is rather pronounced, 
as shown in Fig. 6, it is far from sufficient to form the Z=114 shell gap. It is necessary to point out 
that the isotopic effect is prominent for the two-proton gap. For instance, as observed in Fig. 6, 
the sensitivity of the to the isoscalar-isovector coupling differs clearly for the N=172 isotones 
with Z=114 and 116. Indeed, the pronounced isotopic effect exists for the Z=120 shell closure. 
52p = 3.5 MeV in 292120 reduces to 2.7 MeV in ^^'^12Q. Further, the peak at Z=120 disappears 
totally in ^^^120. As far as the double shell closure in spherical SHN is concerned, ^^^120 turns 
out to be the most possible candidate with various parameter sets in the NL3 calculations. This 
is consistent with the prediction in Refs. [16, 17, 20, 32]. 
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FIG. 7: Single-particle energies in 120 with various parameter sets in the FSUGold calculations. 

Now, we discuss the results with the FSUGold. In Fig. 7, the single-particle energies for 120 are 
plotted with various parameter sets in the FSUGold calculations. Though the isoscalar-isovector 
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coupling is already included in the FSUGold, its strength is changed in FSUGwlS and FSUGw45 
in order to manifest the importance of this coupling. Compared to results shown in Fig. 2, the role 
of the isoscalar-isovector coupling in the single-particle properties is similar in FSUGold, and here 
the nucleus ^^^120 is also doubly magic. Though the role of the isoscalar-isovector coupling is less 
prominent for the empirical shift between 7r3si/2 and 7rl/i9/2 than that in the NL3, it can favorably 
reduce the relative energy between these two levels. Similarly, we can establish correlations as in 
Fig. 4 for the FSUGold results, while here for brevity we neglect the display of the correlations. 
On the other hand, some distinctions of the single-particle spectrum are given with the FSUGold. 
For instance, an inversion of levels ^1^17/2 and ^2/111/2 is observed in the FSUGold calculations, 
while it does not take place in the NL3 calculations. The similar inversion also occurs between 
the iyliii/2 and i>3p orbitals. Moreover, the shell gap at N=:184 is suppressed as compared to that 
with the NL3. On the contrary, the gap for N=198 in FSUGold well develops, and the isoscalar- 
isovector coupling can further enhance the magnitude of this gap. This implies that N=198 is a 
magic number with the FSUGold. With the moderate increase of the proton number, the N=198 
shell gap remains large. For instance, the size of the N=198 shell gap in ^^^126 is even a little 
larger than that in ^^*120. 
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FIG. 8: The same as in Fig. 5 but for the FSUGold calculations. 
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The N=198 shell closure can consistently be observed using the 62n- Fig- 8 displays the binding 
energies and two-neutron gaps for the Z=120 isotopes. As shown in the lower panel of Fig. 8, the 
two-neutron gap at N=198 is clearly increased by softening the symmetry energy. Besides the large 
N=172 gap, we see that the large N=198 gap emerges. The situation of the shell closure at N=198 
in FSUGold differs from that in NL3. This distinction can be associated with different model 
constructions. For instance, the nonlinear self-interaction of the oj meson is included in FSUGold. 
Moreover, the compression modula of the NL3 and FSUGold differ by about 40 MeV. In general, 
the lower incompressibility in FSUGold allows more nucleons to be accommodated in the potential 
well. For the two-proton gap in FSUGold, it is less sensitive to differences in the symmetry energy 
than that in NL3, while the isotopic effect is similar to that in NL3. The increase of the neutron 
number suppresses the Z=120 proton gap. The two-proton gap is 2.7 and 2.3 MeV with N=172 
and 184, respectively. As compared to the NL3 results, the value of 62p with the FSUGold turns 
out to be smaller. The peak disappears with N=198, similar to that with the NL3. 

To examine the consistency of the calculation, we plot in Fig. 9 the change in nucleon density 
distributions in the FSUGold calculations. It is shown that in ^^^120 the neutron density is 
enhanced in the central region. In presence of this central enhancement, the radial inhomogeneity 
appears to increase the sensitivity of the two-neutron gap to differences in the symmetry energy. 
For ^^^120, it was mentioned in Ref. [32] that the magnitude of the central depression increases 
with the decrease of the compression modulus. We note that the central depression in ^'^^120 
with the FSUGold is not more prominent than that with the NL3 as shown in Fig. 3. Indeed, this 
can be attributed to the non-linear w-meson self-interaction in FSUGold that lowers the potential 
barrier. For ^"^^120, as shown in the middle panel of Fig. 9, there exists an inhomogeneity of the 
modification added by the isoscalar-isovector coupling, and this is consistent with the moderate 
decrease of the two-neutron gap at N=184 with the Av, as seen in the lower panel of Fig. 8. 

Next, we investigate the neutron skins in SHN. As it is known, the neutron skin thickness in 
heavy nuclei such as ^osp]-, sensitive to differences in the symmetry energy. In SHN, the case is 
similar. In Fig. 10, we display neutron skin thicknesses in Z=120 isotopes for various parameter 
sets in the NL3 and FSUGold calculations. In general, the sensitivity of the neutron skin thickness 
to differences in the symmetry energy can be well understood in the following way [61, 66, 67]. 
The pressure of pure neutron matter at saturation density is equal to the symmetry pressure. As 
the symmetry pressure decreases with the inclusion of the isoscalar-isovector coupling, the neutron 
skin thickness reduces in neutron-rich SHN. On the other hand, two interesting features in Fig. 10 
are also observed. First, the neutron skin thicknesses in various SHN can roughly be separated into 
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FIG. 9: Nucleon density distributions in the Z=120 isotopes with various parameter sets in the FSUGold 

calculations. 



two reaches according to the slope. The reach with the larger slope is associated with the large 
neutron gap at N=172. With the addition of neutrons above this gap the neutron skin thickness 
thus increases clearly. As the occupation surpasses the much smaller gap at N=184, the nuclear 
attraction from the interior gap can still appreciably restrain the extension of neutrons. This 
leads to a smaller slope at the larger isospin asymmetry. Second, we observe that the difference 
of the neutron skin thickness increases moderately at large isospin asymmetries. Because of the 
looser binding with the increase of the isospin asymmetry, it gives rise to an enhanced sensitivity 
to differences in the symmetry energy. Indeed, the correlation between the neutron skin and 
symmetry pressure in SHN deviates from a simple linear relation for heavy nuclei due to the more 
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FIG. 10: Neutron skin thicknesses for SHN with various parameter sets in the NL3 (upper panel) and 
FSUGold (lower pannel) calculations. 

extended neutron distribution given by a looser binding compared to that of heavy nuclei. 

It is necessary to mention that the uncertainty of the neutron skin thickness is predominately 
from the changes in the neutron radius. Fig. 11 displays the charge radii of the Z—120 isotopes 
for various parameter sets in the NL3 and FSUGold calculations. It is shown that the change in 
charge (or, proton) radii is much smaller than the corresponding neutron skin thickness that is 
shown in Fig. 10. Moreover, for a considerable domain, it is seen that the charge radius of SHN 
can be well approximated as a constant. This is favorable, as one expects that the inclusion of the 
isoscalar-isovector coupling shouldn't compromise the success of models in reproducing a variety of 
ground-state properties. Recently, the precision measurement of the neutron radius of ^°®Pb has 
been proposed at the Jefferson laboratory via the parity-violating electron scattering on neutrons 
in ^°*Pb [68]. The measurement of the neutron radius that promises 1% accuracy will impose a 
strict constraint on the density dependence of the symmetry energy. Correspondingly, this also 
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FIG. 11: The same as in Fig. 10 but for charge radii. 

provides a significant constraint on the properties of SHN through the relationship that can be 
estabhshed for properties between SHN and ^°^Pb. 

At last, we mention a few weak points in this work. First, with the isoscalar-isovector coupling 
included or changed in the best-fit models, we fitted the corresponding parameters without using 
the best-fit procedure. Considering the important effect of the isoscalar-isovector coupling on the 
empirical shift, its inclusion seems necessary in the construction of the best-fit models for the study 
of SHN in future. Second, we note that the effect of the isoscalar-isovector coupling on the empirical 
shift of neutron levels is not as satisfactory as that of proton levels, though the prediction on the 
N=172 shell closure in SHN is not much affected by this coupling. For a more detailed investigation 
of level shifts in future, it may be favorable to consider the coupling with the surface vibration 
modes [69, 70] and the influence of the relatively small Lorentz mass of nucleons in RMF models. 
Third, in this work we have used the simple BCS theory with phenomenological pairing gaps. For 
this point, we address it in some detail. It would be interesting to treat the pairing interaction in a 
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dynamical way such as in the relativistic Hartree-Boguhubov (RHB) theory (for reviews, e.g., see 
Refs. [33, 53, 71]) for the open-shell SHN that are usually deformed, e.g., see Refs. [35-37], though 
the effect of pairing interactions on the shell closure is small. To our knowledge, a RHB model in 
the coordinate space is still not available for deformed nuclei, especially SHN [71], and in the RHB 
model the calculation for deformed nuclei is usually performed with the harmonic oscillator basis 
expansion. Due to the numerical complication, the RHB calculation in the deformed framework is 
still limited, see Refs. [34, 72] and references therein. In fact, most works treat the nucleon pairings 
in open-shell and deformed SHN using the BCS theory. It is interesting to note that the RMF 
model plus the BCS pairing works quite well for the open-shell nuclei except drip- line nuclei [73, 74], 
since for most isotopes the nucleon occupation number in the continuum is just moderate. For 
most open-shell SHN, wc find that the occupation in the continuum is also comparatively small. 
In this sense, the results obtained for most open-shell SHN with the RMF model plus the BCS 
pairing are comparable to those obtained with the RHB model. Recently, a separable pairing force 
was proposed for the RHB model with considerable reduction of the computing time [72], and it 
would be hopefully developed to study properties of open-shell SHN. 

IV. SUMMARY 

In summary, we have investigated the dependence of the ground-state properties of spherical 
SHN on the density dependence of the symmetry energy within RMF models. The various density 
dependences of the symmetry energy are simulated by changing the strength of the isoscalar- 
isovector coupling in RMF models (the NL3 and FSUGold). It is found that the isoscalar-isovector 
coupling produces an important effect on the empirical shift of spherical orbitals in SHN. Especially, 
the empirical shift between the 7r3si/2 and TTlhg/2 in NL3 is nicely reproduced by including the 
isoscalar-isovector coupling. This provides a favorable support for the Z=120 shell closure. In 
addition, the isoscalar-isovector coupling can produce a small but favorable effect on the N=172 
shell closure. With both models NL3 and FSUGold that differ in the compression modulus, the 
double shell closure is predicted in ^^jl^O. The shell closure is also investigated in more Z=120 
isotopes. We have discussed the association of the central depression or enhancement with the 
effect of the softening of the symmetry energy in SHN. Compared to the moderate modification to 
single-particle energies and shell gaps, significant reduction of the neutron skin thickness in SHN 
is expectedly obtained by softening the symmetry energy. Moreover, the proton radius is little 
changed, similar to the situation in ^''^Pb. 
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